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ABSTRACT 

We describe the statistical properties of light rays propagating though a random sea 
of gravity waves and compare with the case for scalar metric perturbations from density 
inhomogeneities. For scalar fluctuations the deflection angle grows as the square-root 
of the path length D in the manner of a random walk, and the rms displacement of 
a ray from the unperturbed trajectory grows as D^/"^. For gravity waves the situation 
is very different. The mean square deflection angle remains finite and is dominated 
by the effect of the metric fluctuations at the ends of the ray, and the mean square 
displacement grows only as the logarithm of the path length. In terms of power spectra, 
the displacement for scalar perturbations has P{k) ocl/k'^ while for gravity waves the 
trajectories of photons have P{k) oc 1/k which is a scale-invariant or 'flicker-noise' 
process, and departures from rectilinear motion are suppressed, relative to the scalar 
case, by a factor ~ (X/D)^^^ where A is the characteristic scale of the metric fluctuations 
and D is the path length. This result casts doubt on the viability of some recent 
proposals for detecting or constraining the gravity wave background by astronomical 
measurements. 

Subject headings: Cosmology - gravitational lensing - gravitational radiation - astrometry 
- galaxy clustering 



1 



1. Introduction 

The deflection of light by scalar metric pertur- 
bations generated by density inhomogcneity is well 
understood: light propagates through an inhonioge- 
neous universe much as it would through a block 
of glass with inhomogeneous refractive index n = 
1 — 2(p — 1 + h, and this provides important observa- 
tional constraints on the distribution of dark matter 
and associated metric fluctuations /i on a very wide 
range of scales. In /i-lensing (angular separations 9 ~ 
10~^arcsec (Paczynski, 1986)) one is probing metric 
fluctuations h ~ ~ 10~^^, and with sufflciently 
compact and distant sources this can be pushed to 
much weaker levels: For femto-lensing of cosmologi- 
cally distant gamma-ray bursters by comet mass ob- 
jects (Gould, 1992) the sensitivity is /i ~ 10~^° local- 
ized in a region of only ~ 10^ cm in size. At the other 
extreme we have weak lensing by large-scale struc- 
ture (Blandford et al, 1991, Miralda-Escude 1991, 
Kaiser, 1992), where the shear 7 and fractional am- 
plification SA are on the order of ~ D / drV^$ where 
D is the path length, and give rms shear and amplifi- 
cation 7, <5A ^ /t(_D/A)"^/^ so the dimensionless shear 
is larger than the dimensionless metric fluctuation by 
a large factor (being the combination of a 'lever-arm' 
term D/X and a 'root-N' factor due to iV ~ D/X in- 
dependent structures along the line of sight adding in 
quadrature). 

It is natural to ask whether one can place similar 

constraints on 'tensor' metric fluctuations — i.e. grav- 
itational radiation — particularly since the most pop- 
ular theories for structure formation such as inflation 
and topological defects make fairly specific predic- 
tions for the gravity wave background; (Abbott and 
Wise, 1984; Battye, et al, 1996, see Allen, 1996 for 
a recent review). For large-scale structure, the cur- 
rent amplitude of tensor modes is very small com- 
pared to the scalar fluctuations, but on small-scales 
the tensor fluctuations are expected to dominate. To 
take an example, popular inflationary models predict 
f^GW ~ H'^lo'^/Hq ~ 10"-^^ on small scales; those 
which re-entered the horizon in the radiation era, cor- 
responding to rms metric fluctuations h ~ 10~^Ho/u), 
where lo is the frequency of the waves and Hq is the 
Hubble parameter. Naively applying the weak lensing 
formula above (with D ^ l/-ff) one finds that while 
the metric fluctuations decrease with increasing spa- 
tial frequency, the amplification and shear increase as 
v^, and reach unity at w ~ lO^^H ~ 10~^Hz. Scalar 



metric fluctuations of this amplitude would certainly 
yield very interesting effects since 6A ^ 1 signals the 
onset of multipath propagation. This would cause 
scintillation and observable time variation of distant 
sources fiuxes (provided the sources have angular size 
6 ^ 10~^^ ~ 10~''arcsec), and there would also be po- 
tentially interesting interference effects at radio wave- 
lengths. However, it is not at all clear that tensor 
metric perturbations deflect light in the same way as 
scalar perturbations. For one thing the perturbations 
are moving, and also the waves are transverse as com- 
pared to scalar modes which cause a longitudinal de- 
flection. As we shall see, both factors have a profound 
influence on the character of photon trajectories. 

This subject has a long but somewhat confusing 
history. Some analyses have assumed that gravity 
waves behave much like scalar perturbations, but with 
rapid time variation: Bergman (1971) discussed scin- 
tillation from nndti-path propagation. It has been 
suggested that gravity waves could cause a time- 
varying lateral displacement of the network of ^- 
Icnsing induced caustics in lenses such as 0957-1-561 
(though Linder (1987; cited in Linder 1988) has sug- 
gested that the effect is rather weak, corresponding 
to an observer-caustic velocity of order ~ h) and 
Fakir, 1994b has discussed the analogous effect for 
caustics produced by ISM refractive index inhomo- 
geneity, see below. Marlcau and Starkman (1996) 
have discussed image broadening by gravity wave in- 
duced 'seeing'. Linder (1988) has calculated the ap- 
parent clustering of galaxies induced by amplification 
by gravity waves. He finds an angular correlation 
function ^0(9) ~ D'^w'^h?, which is what one would 
expect if the rms amplification were on the order of 
5 A ~ hD/X. This is somewhat different from the 
case for scalar metric fluctuations where one flnds 
5 A ~ /i(Z)/A)^/^, so according to Linder the rms 
amplification grows with increasing path length, but 
more slowly than would be the case for scalar pertur- 
bations. However, Zipoy and Bertotti (1968) have ar- 
gued that to first order in the gravity wave amplitude 
there is no amplification, which seems to contradict 
this. 

Braginsky et aL(1990) have calculated the phase 
shift of an initially planar EM wave propagating 
through a stochastic GW background. They find 
(their equation 4.1) a mean square phase shift for 
receivers separated by baseline L and receiving EM 
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waves of frequency uj^^ L 

{{5<i>f)r.,J^^L\hl^ + hl^,) (1) 

This would imply a rms deflection of the light rays 

59 = 54>/{ujryL) on the order of the rnis metric fluc- 
tuation h, but with contributions from the observer 
and source adding in quadrature. A key point here is 
that the phase shift or ray deflection does not increase 
with path length, in stark contrast with the result for 
scalar perturbations where the transverse photon mo- 
mentum performs a random walk and 56 ^ hyJD/X. 
The appearance of a term involving hem raises the 
interesting possibility that one might see relative mo- 
tions of sources which appear close together on the 
sky, and Bar-Kana (1996) in an independent analy- 
sis has obtained a similar result: 50 ~ u)h, and has 
suggested that with VLBI limits on relative motions 
of distant sources this could provide useful limits on 

The works cited above invoke a statistically homo- 
geneous and isotropic stochastic background of grav- 
ity waves. There have also been analyses of the deflec- 
tion of rays passing by a point source of gravity waves 
such as a binary star system. Fakir 1994b has consid- 
ered the case where light rays pass close to a binary 
GW source (impact parameter h ^ Aqw) and finds 
a deflection of order h{b). Durrer (1994) has made a 
similar analysis and finds 56 h{b) in general. 

In an attempt to clarify the situation we will ana- 
lyze the deflection of light rays propagating through 
a stochastic background of gravity waves by means of 
the geodesic equation. We recover the result of Bra- 
ginsky et al., but we show that the term involving the 
metric fluctuation at the observer does not give rise 
to observable image motions. We find that relative 
motions of distant sources are 56 ~ {h/D)y/ln{D/X), 
much smaller than found by Bar-Kana, and conse- 
quently VLBI observations of relative proper motion 
cannot usefully constrain ^Iqw- We find that the 
dominant source of proper motions is the gravity wave 
at the loc;ation of the observer (which gives rise to a 
characteristic distortion of the sky as discussed by 
Fakir (1992) and Pyne et al.{1996)). Aside from this 
we find that the light rays are remarkably straight 
— large-scale deviations from rectilinear motion are 
smaller than the case for scalar perturbations by a fac- 
tor ^ (A/D)^/^ and take the form of a 1-dimensional 
'flicker noise' process (see Press, 1978 for a review of 
flicker noise processes). For gravity waves, the dis- 
placement of a ray grows only logarithmically with 



path length rather than as D"^/^ for scalar perturba- 
tions. We show that stochastic gravity waves do not 
appreciably displace microlensing caustics (in general 
agreement with Linder (1987)), neither do they cause 
galaxy clustering (counter to Linder (1988)), nor do 
they induce observable distortion or rotation of dis- 
tant objects. 

2. Light Deflection by a Planar Perturbation. 

We now compute the deflection of light by a sin- 
gle planar (though not necessarily sinusoidal) metric 
fluctuation. As we will work in the linear approxima- 
tion, the deflection for the general case can be made 
by superposition of planar disturbances. 

We consider geodesies in a linearly perturbed Min- 
kowski spacetime, g^^, = rj^i, + h^j^, with the back- 
ground metric 

= rjf.^dr^dr" = -dt^ + dx^ + dy^ + dz^ (2) 

The generalisation to an expanding FRW cosmology 
is straightforward, but our main point (the qualita- 
tive difference between scalar and tensor fluctuations) 
is adequately illustrated by the simpler flat-space cal- 
culation. 

Consider a photon with 4-momentum p'^ = + W 
where n" is the constant unperturbed vector n'^ = 
(l,n) with n • n = 1 and h" is the first order per- 
turbation due to h^i,. We can obtain the linearized 
geodesic equation by varying the action obtained from 
the Lagrangian 

L = = (j?^, + h^.)p^p''. (3) 

(or by variation of the path length 5 J dt^/L = 0) to 
obtain 

6« = -rj'^^in^hp^ - ^nVdffh^,) (4) 

We take the time coordinate as the affine parameter, 
so the "dot" operator is (d/dt) = (n^(?p). Here we 
are interested in the displacement d of a light ray in 
the plane perpendicular to the unperturbed ray, and 
its time derivative d = b — (n • b)n = b^ which is the 
first order perturbation to direction of the light ray. 

For scalar perturbations we can write /ia/3 = h5a0 
with h = —2$, where $ is the Newtonian gravita- 
tional potential, and we find for the spatial compo- 
nents of the photon momentum bi = —hni +dih. For 
a static planar perturbation with wave direction m 
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we have h = /(m ■ r) so h — nimif , dih ~ m;/', and 
therefore 

d^^h (5) 
n • m 

where mj^ = m — (n • m)n. 

To obtain the corresponding equation for gravity 
waves it is convenient to temporarily adopt a co- 
ordinate frame such that m = (0,0,1) and n = 
(a/I — /i^, 0, /i). The non-zero components for a trans- 
verse traceless gravitational wave are then 



^xx — ^yy — ^+(^ ^) 
^xy — hyx — ^) 



(6) 



and replacing partial derivatives with total time deriva- 
tives much as above (but now with ft, = /(t — m • r) so 
dih = -mif, h = (1-//)/', so dih = -m/(l-^)"^/i) 
we find 

by^-^T^h^ (7) 
k = -^{l + fi)h+ 

The basis vectors in this frame are given by ic = (n — 
(m • n)m)/A/l — /i^, y = n x m/i/l — /i^. z = m, 
so we can write (^ in vector form as h = b^x + . . ., 
and projecting b onto the plane perpendicular to n 
we have 



n X mhy 



(8) 



This is a vector equation and is therefore valid in an 
arbitrary coordinate frame, though it should be kept 
in mind that the bases for the decomposition into po- 
larization states are still defined by n, m. 

The component drives a deflection in the plane 
defined by n, m, and hx drives a deflection perpen- 
dicular to both n, m. Note also that the lengths 
of the vector coefficients of are not equal: 

i|(l — n • m)mj_| = i(l — /x) i/l — while |n x m| = 



y/l — fi^, so for modes with m close to n the coupling 
to the polarization is relatively suppressed. 

A word is in order on the physical meaning of the 
coordinate system we have implicitly adopted here. 
Since /iqq = 0, it follows from (^), (^ that observers 
with f" — (1, 0, 0, 0) initially will remain at rest in this 
coordinate system. Thus equation (^) describes the 
deflections of rays relative to a coordinate grid that 
can be realized physically by a dust of freely falling 
test particles. This is convenient since, for the most 
part, the sources and observers we will consider below 



are in free fall, and therefore coincide with our coor- 
dinate frame aside from some uniform unaccelerated 
peculiar motion. 

To summarize, for planar perturbations, the geodesic 
equation for the transverse deflection of the ray is 



d = T(m, n)h scalar 
d = Ta(m, n)ha tensor 



(9) 



where, for tensor perturbations, there is implied sum- 
mation over the index a = +, x , and where 



T = m^/(n- m) 
and 

T+ = -i(l-n-m)mj 
Tv = — m X n 



(10) 



Thus, for both scalar and tensor perturbations, the 
transverse acceleration can be expressed as the to- 
tal time derivative of the metric times a simple vec- 
tor valued function of the photon direction n and the 
wave direction m. We can trivially integrate d to find 
e.g. the change in the ray direction as the difference 
in the metric perturbation between the end points of 
the ray. To obtain the deflection for a general pertur- 
bation we need to sum over plane wave components. 
This we will do below, but first we note the qual- 
itative difference between the functions T, Ta', the 
former has a pole at = n • m = whereas Ta has 
none. Thus for scalar perturbations modes with ~ 
(which are nearly transverse to the line of sight) will 
have a special significance. This is easy to understand 
physically: For these nearly transverse scalar modes 
the photon stays in phase with the perturbation over 
extended distances and receives a coherent accelera- 
tion over an extended time — these are the 'resonant 
modes' in the language of Braginsky et al. — and the 
ray acquires a deflection proportional ioh/ fi. For sta- 
tistically isotropic random scalar fluctuations, only a 
small fraction of the power is in the nearly resonant 
modes, yet when we compute the variance in the de- 
flection we find that these dominate and give a 'ran- 
dom walk' for the rms deflection c?rms ^ ^JD/Xh ^ h. 

Tensor perturbations are qualitatively quite differ- 
ent. Since the wave crests are moving at c the = 
modes are no longer resonant. The resonant modes 
are now /i = 1, where the photon surfs along with the 
disturbance, but these produce no deflection due to 
the symmetry of a gravity wave and the deflection for 
modes close to /i = 1 is also small. We can immedi- 
ately see that for gravity waves < 1 so the net rms 
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deflection can be no larger than the rms metric fluc- 
tuation, and this immediately rules out the possibility 
that the deflection grows as a random walk. 

3. Deflection by Random Perturbations 

We now explore the statistical properties of rays 
propagating through a general random background 
of metric fluctuations. We flrst give a rigorous cal- 
culation of the 2-point function and power spectrum 
for the transverse deflections of a ray. We then give a 
heuristic derivation. To start, without any loss of gen- 
erality, we decompose an arbitrary metric fluctuation 
into plane sinusoidal waves: 



+ c.c. 



scalar 



hair, t)^y ^/ijk)e'('"*~'^-"-) + c.c. tensor 

(11) 

If we now choose spatial coordinates such that the 
unperturbed ray lies along the z-axis; n — (0, 0, 1) 
and set m = k = cos (p, ^/l^^Jt? simp, n) 

then we can write the geodesic equation as 



d(z) = H(z) 



(12) 



where 



H(2) = y j§:;^h{k)T{k)e'f'''' + c.c. scalar 

H(^) = I / (0/ia(k)T„(k)e»(i-'^)'=^ + c.c. tensor 

(13) 

and where the vectors T, (which lie in the x, y 
plane) have components 



4(l-Ai)yW 



cos if 

sin if 

cos If 

sin ip 



(14) 



smiy9 
- cos ip 



For a static field h{r) we can simply take h{k) to be 
the fourier transform /i(k) — /i(k) = J d'^rh{r)e~^^''^ , 
since this implies /i*(k) = /i(— k) and (^l]) is then 
just the inverse fourier transform. If the field h(r) is 
statistically homogeneous and isotropic, so the auto- 
correlation function for the metric fluctuations = 
{h(r')h{r' + r)) is independent of r' and depends only 
on the modulus of r, then we find 

{h{k)h*{k')) = (27r)3^(k - k')P(fc) (15) 

where P{k) — J d'^r£^fi{r)e~'^^''' is the power spec- 
trum. At the two-point level then, fourier modes at 



distinct k's are uncorrelated. However, when comput- 
ing the mean square values for observables using (11) 
one must allow for the symmetry between h(k) and 
h*{~k), which imply {h{k)h{k')) = {h* {k)h* (k')) = 
{2TrfS{k + k')P{k). 

For the dynamic radiation field ha{r,t) things are 
slightly different. Here the amplitudes for the modes 
in ( |rT| ) are related to the fourier transform of the field 
by h{k) = l{h*(k) + h'*(k)/ik) and h*{k) = ^ih{k)- 
h'{k)/ik), where we have suppressed the polarization 
subscript, and where hik) is the fourier transform of 
/i at i = 0, and h' is the fourier transform of dh/dt, 
also at < = 0, and we now find 

{ha{k)hl{k')) = (27r)3,5(k - k')5abP{k) (16) 

but, in contrast to the static case, there is now no 
correlation between h{k) and h{—k). In both static 
and dynamic cases the total metric field variance is 
(h^) ^ Jdk PP{k)/{2TT^). The Sat in (|l6|) is a conse- 
quence of the assumed statistical isotropy of the pro- 
cess. 

Equation ( [l6| ) clearly obtains for gravity waves 
from inflation where the perturbations were initially 
zero-point oscillations of the metric, so distinct po- 
larization states and wave vectors are completely un- 
correlated, and h{r) is a gaussian random field, but 
they are also valid for a background of radiation from 
coalescing binaries or from decay of cosmic strings 
(Battye et ai, 1996), provided only that the sources 
are randomly distributed in space and in orientation. 
We should stress that ( |l6|) is valid even if the metric 
fluctuations are highly non-gaussian, as, for instance, 
in the case of a background consisting of pulses which 
rarely overlap, though for these types of backgrounds 
there may be non-trivial higher order correlations. 
We shall restrict ourselves to computing variances of 
observables, for which (^^ provides a full description. 

To compute e.g. proper motion of sources etc. it 
is necessary to supply appropriate boundary condi- 
tions and then solve (|l^) to obtain the observable of 
interest (whose variance can then be computed as an 
integral involving P{k)). This we will do presently 
for various interesting observable quantities. How- 
ever, the properties of the trajectories of photons in 
transit are fully determined by the statistical proper- 
ties of the field H(z) which is just some statistically 
homogeneous 1-dimensional random process. To ex- 
plore this we now compute the spatial 2-point func- 
tion for the H field and its analogue in fourier space, 
the power spectrum. From (13), (M), (fq) we find 
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{R„MHn{z)) = \5„,^S,h{z) with 

Ch(z) = (H(0).H(z)) 



= \!4^P{k){Tl+Tl) 





+1 
-1 



4\ i{l-ti)kz 



(17) 

from which we can readily obtain the power-spectrum 
for the H-ficld 

PH{k) = / dz^H{z)e'''' 

OC 

^ / dk'k'Pik') 



1 

X / dn{5 - 2^ - 4^2 ^ _ ^/)^(^ _ (1 _ fc/A:')) 
-1 

oo 

^ _|_ J dfc'F(fc')(8 - 4/c//c' + 2(A:/A:')^ - (fc/fc')^) 

fe/2 

(18) 

For scalar perturbations we find from (|l|). 



^h(2) = 



47r2 



dkk^P{k) / - ^2)^-^6*^'=^ 



(19) 

Interestingly, ^^(2) is not well defined for scalar 
perturbations as the /i-intcgral diverges at /i ^ 0. 
However, the 'structure function' Sh{z) = (|H(0) — 
H(z)|2) = 2{£,h{Q) - iniz)) is weh defined, as is the 
power spectrum 



P^(fc) = ^ ] dk'k'P{k') J dn{l - /i2)^-2^(^ - k/k') 

-1 



= J dk'k'^p{k'){i - 

(20) 

though this is divergent in the limit fc ^ 0. 

The power spectra (|l^, |o|) are quite revealing. The 
field H(z) is, aside from a constant determined by 
boundary conditions, equal to the direction of the ray, 
so Pj^{k) — Pf/(fc) and the power spectrum for the 
displacement is just Pd{k) — PH{k)/k^. Similarly, 
the power spectrum for the transverse acceleration is 
Pj{k) = k^Pnik). We are primarily interested in the 
deflections and displacements of the beam over long- 
path lengths D ^ Xh, where is the characteristic 
wavelength of the metric fluctuations. These are de- 
termined by the long- wavelength behavior of -Pj(fc) at 



k <^kh, in which case we can drop the terms involv- 
ing k within the integrals and also relace the lower 
limit on the /c^-integration by zero and find the lead- 
ing order term for the displacement power spectrum 



Pd{k) = 



-i J <^k'^P{k') 

dk' 



{h'^)kfi/k'^ scalar 
J ^P{k')^ {h^)/{klk) tensor 



(21) 

Thus the effect of both types of metric fluctuations 
is to produce a universal scaling law for the power- 
spectrum of displacements at fc <C fc/i, but the form 
of the scaling is quite different: Scalar modes give 
a displacement which is the double integral of a 
'white-noise' process: Pd{k) cx fc^**, and the mean- 
square deviation, averaged over a scale D is (d^) d 
{kPd)k=i/D OC D^. Tensor modes, on the other 
hand, generate a '1//' or 'flicker noise' spectrum 
Pd{k) cx k^^ with scale-invariant deviation {d"^) ~ 
{kPd)k=i/D ~ ^hi^^) ^ P*^- For a given scale and 
amplitude of metric fluctuations the large-scale de- 
viations from rectilinear motion are suppressed by a 
factor '--^ (A/Z?)'^/2 foj- tensor modes as compared to 
scalar modes. 

Another interesting feature emerges if we consider 
the power spectrum for the acceleration at zero fre- 
quency. In the scalar case P^{k) tends to a constant 
value at fc ^ whereas for tensor fluctuations P^jik) 
being proportional to k^ vanishes in this limit. How- 
ever, the power at zero frequency is the integral of the 
corresponding 2-point function: -Pj(O) = / dz^^{z). 
What this is telling us is that for scalar perturbations 
the sequence of accelerations suffered by a photon 
can be legitimately modeled as a series of uncorre- 
lated kicks (as this gives the correct random walk for 
the photon direction) whereas for tensor fluctuations 
the accelerations are strongly anticorrelated, and this 
again helps us understand why the net deflections are 
so much smaller for gravity waves. 

Equation ( pl| ) give the universal asymptotic scaling 
laws for the low-frequency power spectrum of photon 
trajectories propagating through metric fluctuations 
of relatively higher spatial frequency. They are the 
main result of this paper and show immediately that 
tensor metric fluctuations are extremely inefficient at 
deflecting light. 

We can give a somewhat simpler derivation of these 
scaling laws. The mean-square systematic deflection 
of a ray over path length Z), is proportional to the 
power in spatial modes which project to low spatial 
frequency ~ 1/D along the line of sight, times the 
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square of the amplitude factor T. For scalar per- 
turbations, if we compute (d^)D ^ {kPj^{k))k^ii £, 
we are only sensitive to modes which have /i 
l/{khD), where kh is the spatial frequency of the 
metric fluctuation. These are a fraction ^ ^ oi 
all modes, and have an amplitude factor T ~ l/jj, 
so {d^)D - {h^){khD)'^T^ - khD{h^) and hence 
{<P)d ~ D^{d^)D ~ khD^{h^). For tensor fluctua- 
tions, the modes which project to line-of-sight wave- 
vector k ^ 1/D are those with (1 — /i) ~ l/fccw-D. 
These again are a fraction ~ l/fccw^ of all modes, 
but they have T+ ~ (kowD)-^/^, ~ (fccw-D)"^/^ 
and consequently, for k^y^D ^ 1, the x -polarization 
components dominate, {A^)d ^ {h'^){kGwD)~^T^ ~ 
{h'^)/{kGwDf and therefore (d^)^ ~ ^'^(d^)^ - 
k^^{h?) which is independent of scale D. 

We will now work through some calculations of 
what are, in principle, observable quantities. From 
the order-of-magnitude estimates above this is largely 
an academic exercise, but will reveal some subtleties 
in handling boundary conditions. 

4. Applications 

We will now apply the results obtained above to 
compute the deflection of a coUimated beam (§4.1.); 
prope r mo tions of sources (§ [4.2j ); deflection of caus- 
tics (§ [4.3j ); and amplification, distortion and rotation 
of images of distant galaxies (§4.4.) 



4.1. Deflection of a Collimated Beam 

Let us calculate the change in the direction d for 
a ray which propagates from zq to zi = zq. Clearly, 
from (|l|) we have 



so 



M = d(zi) - d(zo) = H(zi) - H(zo) 



(|<5dp)^2(^H(0)-eH(i?)) 



(22) 



1 

Terr' 



/ dkk^P(k) 




+1 

X / d^l{f> 



(23) 

where D = zi — zg is the path length. For kD ^ 1 
we can neglect the oscillatory term in the /.t-integral 
— i.e. ^h{D) ^ C-f/(0) — and so the integral is ele- 
mentary and we have 



15^ ' 



This is in accord with the order of magnitude estimate 
of Braginsky et ai. Clearly, the directional change 
does not grow with increasing path length. Compare 
this with the analogous situation for scalar fluctua- 
tions where we have 



(IMP) = ^J dkk'P{k) J d^il-i^(l_e^Mfei.) 

-1 

(25) 

for fixed k and for D ^ 1/fc the fi integration is now 
dominated by /i ~ l/{kD) (these are the 'resonant 
modes' where the photon sees a phase change over 
the path length on the order of unity) and we have 



(IMP) dkk^P{k) J dyy-Hl - e^y) 



(26) 



- Dk^{h^) 

where we have defined the characteristic wave-number 
_ / dkk^P{k) 



jdkk'^P(k) 



(27) 



(24) 



This is very different from ( p4| ) and the variance in the 
deflection angle now grows linearly with path length 
as expected for a random walk. 

Note that the ratio of (||) to (|2|) is not - D^/\l, 
as one might have anticipated from the discussion of 
the previous section. This is because ( |2^ ) is domi- 
nated by 'surface terms', rather than low-frequency 
power with kh ^ 1/D. 

4.2. Proper Motions 

In the preceding section we found that the mean 
square change in direction of a ray propagating through 
random gravity waves is essentially the sum of the 
mean square metric fluctuations at the end points. 
Does this mean that two sources close together on the 
sky but at different distances would show a relative 
proper motion on the sky with amplitude S9 huj? If 
so, then with VLBI observations one could constrain 
^GW as described by Bar-Kana, 1996. Unfortunately 
this is not the case. The foregoing correctly gives the 
final direction of a collimated beam which initially has 
d = 0. Physically, one could realize this with a pair 
of neighboring freely falling observers, one of whom 
shines a flashlight at the other who is holding a plate 
with a pin-hole in the centre which selects a particular 
beam. This is rather like a lighthouse, where a beam 
suffers a deflection at the source, but where distant 
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observers do not see a moving source, rather they see 
a flashing source (and whenever they actually see the 
source, they see it at the same location). 

Clearly, to compute how real astrophysical sources 
— which are generally at best poorly coUimated — 
would move we need to 'solve the lens equation' and 
find the final direction for a ray which leaves the 
source at zq and actually arrives at the location of 
the observer d(zi) — 0. This is readily done. Inte- 
grating (n2|) once more we find 



dzB.{z) (28) 



d(z) = (z - zo) d(zo) - H(zo) 



so the beam which reaches a detector at d(zi) — 
must set out with 



d(zo) - H(zo) 



1 



zi ~ zo 



dzH(z) (29) 



and consequently will have a direction at the observer 



d(zi) =H(zi) 



1 



zi - Zo 



dzH(z) 



(30) 



from which we see that the effect of the metric fluctu- 
ation at the source enter only as a minor contribution 
to the integral of the metric fluctuations along the 
line of sight. 

Clearly, for a random fluctuating H(z) and dis- 
tant sources the second term will be sub-dominant. 
However, the first term depends only on the met- 
ric fiuctuations at the point of observation, and is 
a slowly varying function of the angle between the 
wave-direction and the linc-of sight. This term causes 
a characteristic distortion of the sky which is the same 
for all distant {D ^ A) sources. This is the effect dis- 
cussed by Fakir (1992), Pyne et a/.(1996). One has 
to be careful in how to interpret this term as we have 
derived it. The direction here is defined as relative to 
our coordinate system which is tied to freely-falling 
test particles. Thus d here would correctly describe 
the apparent motion of a distant source relative to 
a nearby {D <C A) freely-falling foreground reference 
source, but would not directly describe the change in 
relative angles between distant sources as measured 
with a sextant. To compute the latter we need to 
calculate how our nearby reference source would ap- 
pear to move in the physical reference system. This is 
straightforward, and gives a somewhat different pat- 
tern of displacements (and we recover their equation 
51), but does not change the general result that the 



angular deflection is on the order of h and is slowly 
varying across the sky (and so would not be detectable 
by measuring relative proper motions of neighboring 
sources with VLBI). 

Let us now estimate the size of the second term in 
(^), which depends on the metric fluctuations along 
the line of sight to the source, and which would there- 
fore be expected to give rise to relative motions of 
sources which appear close together on the sky but 
lie at different distances. The variance of the integral 
J c?zH(z) is 




dzH{z] 



sin^ kD/2 



(31) 



Now from ( [T^ ) we see that Pnik) consists of four 
terms of the form fc" j^^^^ dk'k'^~'"P{k') for n = 
1 ... 4, each of which rises as a power law for k ^ /cqw 
but then falls to zero for k 3> fccw as the lower 
limit on the integral takes effect. It is not diffi- 
cult to see that for kcwD 1 the leading contri- 
bution to the variance above comes from the first 
term Ph ~ k J^^^^ dk'P{k') and the integrand in (31) 



J dk'P{k')/kD^ and so there is equal contribution 
to the variance above for each logarithmic interval in 
the range D^^ <C fc ^ fccw (with the sin^ kD factor 
limiting the divergence at low k and the lower limit 
on the J dk'P{k') integral limiting the divergence for 
k ^ fccw) a-iid 




dzH{z) 




Xl^liiiD/XGw){h^) (32) 



which is just what one would anticipate from the gen- 
eral idea that the photon trajectories are a cx 1/fc 
or fiicker noise process. To obtain the mean square 
angular deflection we simply divide (|3^ ) by the source 
distance D^. 

Thus the net effect of the the perturbations along 
the line of sight is to displace the image by an angle 9 
on the order of {X/ D) ^J\n{D / X)h which is very small 
indeed. Of course, as we do not know the unperturbed 
image location this is not directly observable, but its 
time derivative 9 = 89 /dt would be observable as a 
proper motion (or as a relative proper motion for two 
sources which are close together on the sky but lie at 
different distances) . It is straightforward to compute 
(^^), since the time variation simply introduces an 
extra factor fc^ in e.g. (18) so the power spectrum for 



8 



dH/dt is ~ fc / dkfikf-i^P{kh), and wc find 

(^2)1/2 ^ {h/D)^/HD/X^ ^ nl^^HoXaw/D 

(33) 

For gravity waves with periods 1 / uj 10 yr consid- 
ered by Bar-Kana (1996), and for cosmologically dis- 
tant sources the rms proper motion is 56 ^ 10~^^\/riGW 
arcsec/yr which is 9 orders of magnitude smaller 
than obtained by Bar-Kana, and would be rather hard 
to measure. Similarly, one can estimate the broad- 
ening of images due to gravity waves (which must 
have periods less than the integration time tint) to 
he 9 ^ ^^o^fntV^GW which would be on the order 
of 10~^^-\/f^GW arcsec for integration times of a few 
hours. 

We should emphasize that this is only the relative 
motion of sources. The first term in ( |30| ) gives a much 

1/2 

larger motion 9 ~ ujh fiJ^iJo, but as discussed 
above this is a purely local distortion effect and would 
require precise relative astrometry over large angular 
scales, and would not be easily accessible with current 
technology. 

For scalar fluctuations, the variances of H and 
J dzH{z) are both ill defined due to divergent low-fc 
behavior noted above, but the variance of the com- 
bination in (^0|) is however well defined, and {69'^) ~ 
{D / \){h'^) . Note the qualitative difference in the de- 
pendence of the angular displacement on source dis- 
tance; for scalar perturbations the deflection grows 
without limit, while for tensor fluctuations, and aside 
from the local sky distortion offect, the deflection 
tends to zero for very distant sources. 

4.3. Wiggling Caustics 



We found above (4.1.) that a coUimated beam can 



suffer a deflection due to the metric fluctuation at the 
source which leads to a relatively large 'lighthouse' 
effect, but noted that highly coUimated sources were 
rare in Nature. However, a notable exception is the 
caustic network produced by /x-lensing whose sharp 
cusp-like features cause observable time variations of 
source intensities as they sweep past an observer. The 
same is true of caustics arising from refractive index 
inhomogeneity from the ISM. This raises the interest- 
ing possibility that gravity waves, if they can later- 
ally displace the caustics, might give rise to anoma- 
lous time variations of sources which might be recog- 
nizable. Fakir (1994b) has studied the deflection of 
caustics which pass close [b ^ Agw) to a binary GW 
source. Here we will compute the effect of a stochastic 



background. 

We can readily compute the relative motion be- 
tween a caustic network and a freely falling observer. 
Consider a source at zq and a lens at zi , and let both 
of these have d — 0. Now consider a ray which propa- 
gates from the source to the lens. As before, imposing 
the boundary conditions d(zo) = d(zi) = gives the 
direction on arrival at the lens of 



Hi 



1 



dzB. 



(34) 



so, as the photon propagates to z > zi it will have 



d = di+H Hi=H — / dzn 

Dqi 



(35) 



and hence will arrive at Z2 with 



do = / dzd ~ \ dzH — 



Dl2 

Dm 



dzVL 



(36) 



If both Dqi I are ^ Agw and are of similar 

magnitude then d2 is composed of integrals of the 
same form as in the previous section and consequently 
AGw'i\/ln(I?/AGw), which is very small. 



(di)V2 



and the velocity with which the caustic moves across 
the observer plane is u ~ /i'y/ln(Z?/AGw) which, 
aside from the logarithmic factor, agrees with the 
rough estimate of Linder (1987). The exception is 
if the lens is very close to the source, where we can 
get a somewhat larger effect. For Dqi ^ Agw, 

J dzH ~ Z^oiHo and we get a stronger 'lighthouse' 

effect with d2 = Z)i2Ho and caustic-observer velocity 
V = dd2/dt ~ _Di2Ho/Agw ^ £'i2^/Agw- Unfortu- 
nately, however, the cross-section for /i-lensing (of for 
formation of caustics by inhomogeneous ISM electron 
density) vanishes for Dqi — s- 0, and the typical caus- 
tics have i^oi — -D12 and the effect is very small, and 
is certainly much smaller than typical peculiar mo- 
tions of observers and lenses for any plausible value 
of JIgw- 

4.4. Source Amplification, Distortion and Ro- 
tation 

As a final application, we consider the possibil- 
ity of amplification, distortion and rotation of dis- 
tant sources. As discussed in the Introduction, Linder 
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(1988) has suggested that deflection of hght by grav- 
ity waves wih change the number density of sources on 
the sky in an inhomogeneous way and give rise to ap- 
parent galaxy clustering. For static lenses, the change 
in number density is directly related to the amplifi- 
cation of the flux density, since if the lens mapping 
increases the area (and hence flux, as surface bright- 
ness is conserved for static lenses) of a source by some 
factor A, then it will also decrease the expected num- 
ber density of sources by the same factor. This is for 
'standard-candles'. For real sources we have to worry 
about the change in detectability with flux, but in 
general we expect the change in surface density of 
sources to be proportional to the amplification (with 
a constant of proportionality known as the 'amplifi- 
cation bias factor' which can be positive or negative 
depending on the slope of the number-flux relation). 
However, as we also mentioned, Zipoy and Bertotti 
(1968) have argued that there is no amplification due 
to gravity waves, which seems to contradict this. To 
resolve this we will estimate the expected effecting 
using our formalism. We will also discuss the distor- 
tion and rotation of distant sources, which are not 
excluded by Zipoy and Bertotti's argument. 

This calculation differs from the previous applica- 
tions in several respects. First, to calculate the ampli- 
fication, distortion and rotation of sources in a single 
plane we need to solve the geodesic deviation equation 
for a pair of neighboring rays to obtain the mapping 
between angles at the observer and distances on the 
source plane (which we can express in terms of an am- 
plification matrix Sim + ipim where the 'distortion ten- 
sor' is linear in the perturbation amplitude). Sec- 
ond, we then need to average this over source planes 
distributed according to the selection function for the 
galaxies. Also, we need to include the possibility of 
modulation of the surface brightness due to the non- 
stationarity of gravity waves. 

For scalar perturbations, it is well known that 
the mean square amplification is on the order of 
{D/X)^{h'^) (and hence the induced angular correla- 
tion function 'w{9)) grow with depth of survey as D^. 
Linder (1988) found w{9) - {D/Xf{h'^) cx D"^ . How- 
ever, this would be hard to reconcile with the result of 
§ 4.2. where we found that the angular displacement of 
a distant source was composed of two terms: a large 
angular scale distortion of the sky due to the gravity 
wave amplitude here and now, which is on the order 
of ft., and is the same for all distant {D ^ A) sources, 
and a small extra deflection 66 ^ {X/D)^\n{D/\)h 



which depends on the source location, and tends to 
zero for very distant sources. Now the components of 
the distortion tensor ipim are the angular derivatives 
of the deflection angle. Thus, the local distortion ef- 
fect would be expected to produce t/; ~ /i, since the 
coherence scale for this term is on the order of one ra- 
dian. The non-local term due to metric fluctuations 
along the line of sight will cause a fluctuating dis- 
tortion which may be more accessible to e.g. galaxy 
clustering studies. If one imagined that the coher- 
ence scale were 6c ~ X/D i.e. the angle subtended 
by one wavelength at the distance of the source, then 
one would predict distortion 'ip ^ ^Jhi{D / X)h which 
is somewhat larger than the local sky distortion effect 
(though still very small compared to Linder's calcula- 
tion). In fact, as we shall see, even this overestimates 
the true amplitude of the effect. 

Consider two neighboring rays which are nearly 
parallel to the z-axis. Let the separation between the 
rays be 

Ar = Aro + Ad (37) 

where the zeroth order separation obeys Aro — 0; and 
where 

Ad; = Aropdpdi = Arop^p; (38) 

where the latter equality defines the two-dimensional 
transverse tidal field tensor cjjpi . For scalar perturba- 
tions, (f)pi — dpdih, which is a statistically homoge- 
neous random field with fourier decomposition 



where 



(2^ 



k^Rpi(k)h(k)e 



ifikz 



Ri„,(k) = (1 - A*' 



cos^ ip 
COS if sin if 



cos If sm if 
sin^ (f 



(39) 



(40) 



and where we have used dp — > ikp — ■\/l — /i^{cos (p, sin 
For tensor perturbations, we had d = Taha, or from 



d, = 



so the tidal field tensor is 



<^pz = dpdi = -i / ^T;(k)ft(k)(l - M)fcfcpe^(i-^)'=^ -I- 1 
= 5 / x0^fc'-Ra/„^/^a(k)e^(l-^)'=^ + c.c. 

(42) 
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where we have defined 

R+lm = 



cos^ (f COS ip sin ip 
cos (f sin (fi sin^ ip 



cos (/3 sm (p 

■ 2 

sm 95 



-cos^ (p 



-cosp smip 
(43) 

To obtain the shear tensor for a source plane at 
some distance Zg we need to integrate (38) back along 
the null rays to obtain the source plane separation Ar 
for a pair of rays which arrive at the observer with 
angular separation 0. These rays have zeroth order 
separation Aro = 9z, so, with boundary conditions 
Ad ~ Ad = at the observer the solution of ( ^8|) is 
just the particular integral: 



Adi{zs) = Op J dz J dz' z'(j)ip{z') = J dzz{zs-z)(j)ip{z) 
00 

(44) 

and hence obtain the mapping between source plane 
separation Ar = Afq + Ad and angle: 



Ari{Zs) = e.mZs [Sim + i'lm] 



(45) 



where the distortion tensor is a certain integral along 
the line of sight of the tide: 



i'hnizs) ^ J dzz{l~ z/Zs)(t)im{z) (46) 


We are interested in the expectation value for the dis- 
tortion taken over the distribution of distances n{z) 
for the galaxies whose clustering, shapes etc. we are 
measuring, which we can express as 

00 CXD 

^im^ J d,Zsn{zs)ipim{zs) = J dz(f>im{z)g{z) (47) 


with 

00 

g{z) = / dz'n{z'){z - 7? j z') (48) 



for scalar and tensor perturbations respectively, and 
where 

~g(k) = J dzg{z)e'^^ (50) 

We can now compute the variance of the distortion 
matrix elements. First though, it is convenient to 
decompose the 2x2 tensor i/^j^j into a four component 
entity 

Va = Malmi>im. (51) 



with 



Mo = 
Ah = 



1 

1 

1 

1 



Ml = 





I -1 
1 
-1 



(52) 



so Fq is the convergence, Fi, r2 are the shear and F3 
describes rotation of the image. Summing over fourier 
components, we find 



where 



/(|;;f.fc4p(fc)i?„(k)i?Mk)|g(M)P 

(53) 



Ra = {1 — ^^)[l,cos2(y5, sin2i^,0] 
R+c = 5(1 - /i^Kl -m)^[1, cos 2(^, sin 2(^,0] (54) 
i?XQ = (1 ^ ^ /i)[0,sin2</j, cos2(^, 1] 

Performing the (^-integration in (|^ we find that 
the decomposition (|T]) results in a covariance ma- 
trix (FqF^) which is diagonal. For scalar pertur- 
bations F3 vanishes — no image rotation — and 
(F^) = (F2) = 5(1^0) with mean square convergence 



00 ^1 

(r^) = ^ / dkk'Pik) J d^l{l~^lY\~g{pk)f (55) 



To estimate the value of the /i-integral here we need 
to know the asymptotic form of g{k) for large and 
small k. One can see on dimensional grounds from 
(H) ^) that |g(2/)p ~ 1)4 for y ^ To obtain 

the high frequency asymptote we use 



and where we have taken n{z) to be normalized to 
/ dzn{z) = 1 so, on substituting (^^, and (^) in 
(EtI) wc find 



H (l#ffc'^f"(k)/^(k).9(M) + c.c. 

5 / ^fc'^aim(k)/la(k)ff((l - fl)k) + C.C. 

(49) 



g{k) = J dze''^'^ J dz'n{z'){z - z"^ / z') 
2 

00 z' 

= / dz'n{z') J dz{z - z"^ I z')e'^'' 



h J dzniz)il + e"=-) + f / dz^il - e'^^-) 

(56) 
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In general one can consider n(z) to be composed of 
the smoothly varying expectation value n{z) plus a 
fluctuating component due to galaxy clustering. We 
will first consider the smooth term n — n. To obtain 
the leading order behavior of the two fouricr type in- 
tegrals at high k we use the series expansion, obtained 
by repeated integration be parts, 



Z2 JV-1 

J dzfiz)e''^^ = E (-1 

Z2 



/""(z)e"' 
(ife)"+i 

+ ^Idzf^''\z)e^'^^ 



(57) 



which is valid provided the 0th through A^th deriva- 
tives are well defined. Thus, if f{zi) ^ /(^2) then we 
will find a leading order term / oc 1/fc, otherwise we 
have / ~ unless /'(zi) = f'{z2) in which case 

/ ^ and so on. 

Here we have zi = 0, Z2 oo, and we will assume 
that the integrands and all necessary derivatives van- 
ish as z ^ oo, so provided zn{z) and n{z)/z 
is finite as z ^ then the two fourier integrals in 
(56) fall at least as fast as 1/k^ for large k. Thus, 
the dominant term at high frequency is the first term 
g = — and hence |gp ~ for large k. The /u- 
integral in (55) is therefore dominated by modes with 
fj, ^ 1/kD for which (1 — //^) is very close to unity 
and consequently 



(r§) 



Jdkk^P{k) f dy\g{y)\^^ 



D^Jdkk^P{k) - {D/Xf{h^ 



(58) 



which is the usual weak lensing result. 

For the tensor case things are somewhat more 
complicated: Now the leading order contribution to 
\gi{l - /x)fc)|2 in (M) falls as 1/((1 - n)k)'^ for large 



k, but the pair of i?^ and terms contain 6 and 4 
factors of 1 — /i respectively, so now the /^-integration 
is not restricted to small 1 — n and has value ~ 
and hence we find typical covariance matrix elements 
(r^) ~ {h'^), with all four components of Fq hav- 
ing similar mean-square expectation values. However, 
the leading order term in g is also independent of the 
depth of the survey, so what we have calculated here is 
clearly the shear due to the local sky-distortion effect. 
To extract the non-locally generated shear we need 
to subtract from g the distance independent contri- 
bution. We then find asymptotes |g(fc)P ~ l/Zc* for 
k < 1/D and \g{k)\^ - l/{D^k^) for k > 1/D. We 
now find that the dominant contribution to (F^F^) 



comes from the ' x '-polarization and is dominated 
by the nearly resonant modes with 1 — /i ~ 1/kD. 
Now from (^) one can see that Fxo — 0, so the x- 
polarization gives no convergence; Fq = 0, so there 
is no area amplification (in accord with Zipoy and 
Bcrtotti), and (F^) = (F|) = |(F§) with mean-square 
rotation 



in) 



oo oo 

^ I dkkP{k) I dyy^\~giy)\' 



{X/D){h^ 



smaller than the scalar case by 



(59) 
factor 



ri) 



which is 
iX/D)\ 

It may seem a little surprising that {66'^) 
yet the mean-square shear (7^) = (Fj 
{\/D){h^), since, on general grounds, one expects 7 ^ 
59 /6c where Of, is the coherence scale for the angular 
deflections and so it would appear that 9c ~ \/ A/D, 
which is much less than \/ D which is the angle sub- 
tended by the characteristic wavelength or coherence 
scale for the metric fluctuations. The reason for this 
is that the shear is dominated by the nearly resonant 
modes with 1— /i ^ 1 /kD, which run nearly parallel to 
the line of sight, and these have a projected transverse 
wavelength ~ A/ — /i^ VXD, and therefore sub- 
tend an angle 9c ^ ^\/D. 

We have found that the leading order non-local ef- 
fect gives no convergence whatsoever. However, there 
is in fact some flux modulation due to the Sachs 
Wolfe effect, which causes a modulation in the surface 
brightness given by the difference in the metric fluc- 
tuations at the source and observer. One can calcu- 
late the mean square ampliflcation as before by sum- 
ming over plane waves, and again only the nearly reso- 
nant modes contribute and we flnd (Fq)^/^ ^ \J \/ Dh 
which is of the same order as the shear and rotation. 

The shear 7 ~ /i due to the local sky distortion is 
much larger than the spatially inhomogeneous shear 
7 yJXfOh due to the metric fluctuations along the 
line of sight. The gradients of these two contributions 
to the shear are of similar order, however, and high 
frequency fluctuations in the shear are dominated by 
the non-local term. 

We have assumed above that n(z) ~ n(z). Fluc- 
tuations in n{z) can increase the predicted shear and 
amplification somewhat. If we consider the rather ex- 
treme case of a narrow beam which intercepts only 
one cluster (of width bzc) then one finds that (7^) ~ 
(/i^) if bzc < 1/fc and {'-^'^)) / [kSzc) if 5zc > 1/fc, so 



{h''){\/Df\n{D/X) 
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for a cluster which is smaller than the characteristic 
wavelength the shear, amplification etc. can be of or- 
der h. However, this does not make the effect any 
easier to observe since in the situation we have de- 
scribed — which might be a good model for a narrow 
'pencil-beam' survey — the fluctuations in the num- 
ber density from beam to beam would be of order 
unity. 

Our analysis has yielded a very different result from 
Linder's. We have not been able to pinpoint the 
source of the discrepancy, but is would seem to be 
a consequence of inappropriate boundary conditions; 
one can view (^8|) as a second order ODE for ■0im ■ For 
generic initial conditions, one will flnd, in addition to 
the particular integral solution terms which are con- 
stant or grow linearly with source distance, and which 
could plausibly account for Linder's result. However, 
as we have seen, with the correct boundary conditions 
at the observer these terms are absent. 

The weak lensing effects obtained above are very 
small indeed. For the local term, we find 7 etc. on the 
order of ^ ~ \/TIqy/ Hq / lu , and for the non-local 
term, and for distant sources {D ^ ^/Hq) we have 
7 ^ Vf^Gw(-ffo/w)^/^ ^ \/UGwd~^, and therefore for 
waves with 6c ^ 10~^ say we would find an induced 
clustering effect w{9 ~ 9c) ^ 10~^^^2gw which is very 
small for any plausible value for J^gw- 

5. Discussion 

We have calculated the statistical properties of 
photon trajectories propagating through a statisti- 
cally homogeneous and isotropic sea of metric fluc- 
tuations. 

Stimulated by the power of lensing phenomena to 
constrain small amplitude metric perturbations and 
the theoretical expectation that the dominant metric 
fluctuations on small-scales are most likely predomi- 
nantly in the form of gravity waves, we have focussed 
on the effect of metric fluctuations with wavelength 
much less than the path length. 

We have found qualitatively different behavior de- 
pending on whether the metric fluctuations are scalar 
or tensor (gravity wave) fluctuations. In the former 
case, the photon trajectories can be pictured as suf- 
fering a set of uncorrelated kicks which give a random 
walk for the photon momenta (and the displacement, 
which is the integral of the transverse momentum, 
grows as the 3/2 power of path length). For gravity 
waves we have found that this picture of uncorrelated 



kicks is grossly misleading and that the trajectories 
are in fact a scale-invariant 'flicker-' or 'l//-noise' 
process. While a plot of a photon trajectory over 
a short distance of a few gravity wavelengths would 
look qualitatively similar for the two cases, on larger 
scales the coupling for gravity waves is suppressed by 
(A/D)^/^, and the photons can propagate over great 
distances with barely any deflection. 

In addition to this rather general result, we have 
made detailed calculations of a number of potentially 
observable effects: proper motions of distant sources; 
lateral displacement of caustic networks; and weak- 
lensing of distant objects. In all cases the estimated 
size of the effects are too small to be of much practical 
interest to experimentalists, counter to some claims in 
the literature. 

Our analysis has concentrated exclusively on stochas- 
tic backgrounds of gravity waves and so our results 
cannot be directly applied to deflection of light pass- 
ing close to a binary system emitting gravity waves 
(as discussed extensively by Fakir and Durrer). We 
have however, attempted to repeat these calculations. 
For very small impact parameters (6 <, A) we confirm 
the deflection on the order of SO ~ h{b) as found by 
Fakir (1994a,b) (though we note that this effect is on 
the same order as the scalar metric perturbation as- 
sociated with the binary) but that at larger impact 
parameters the effect falls off as a high power of 1/6, 
and we were not able to conflrm Durrer's claim that 
66 ^ h{b) in general. 

To end on a slightly more positive note, we have 
found that the dominant astrometric effect of a grav- 
ity wave background is the local distortion of the 
sky due to the metric fluctuation in the vicinity of 
the observer (see Pyne et al, 1996), which will cause 
proper motions 66 ~ Luh. Barring the possibility of 
a strong background of high frequency waves from 
coalescing binaries, the dominant contribution to this 
effect would most likely come from waves on the order 
of the present horizon size. While the predicted size of 
the motions is impressively small (~ 10~^arcsec/yr), 
and would typically be similar in amplitude to the 
~ 30km/s or so motions due to horizon scale scalar 
metric fluctuations they might in principle be mea- 
surable as a relative motion of distant galaxies. 
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